Stable smectic phase in suspensions of polydisperse colloidal platelets with identical 
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We report the nematic and smectic ordering in a new aqueous suspension of monolayer a- 
Zirconium phosphate platelets possessing a high polydispersity in diameter but uniform thickness. 
We observe an isotropic-nematic transition as the platelet volume fraction increases, followed by 
the formation of a smectic, an elusive phase that has been rarely seen in discotic liquid crystals. 
The smectic phase is characterized by X-ray diffraction, high-resolution transmission electron mi- 
croscopy, and optical microscopy. The phase equilibria in this highly polydisperse suspension is 
rationalized in terms of a theoretical approach based on density-functional theory. 

PACS numbers: 61.30.Eb, 64.70.M-, 81.16.Dn, 82.70.Dd 



Synthetic inorganic platelets are being used extensively 
to investigate phase transitions in discotic liquid crystals 
consisting of colloidal, approximately disc-shaped, parti- 
cles in aqueous suspension. Starting from the completely 
disordered isotropic (I) phase in the dilute regime, most 
discotic suspensions exhibit an orientationally ordered 
nematic (N) phase as the volume fraction of the platelets 
in the sample increases. At higher volume fractions a 
columnar (C) phase, with columns of platelets arranged 
in a two-dimensional lattice, has been observed in some 
materials .either directly from the I phase [l[ or via N 
ordering Q. Computer simulations on monodisperse Q 
and polydisperse [J] hard platelets of high aspect ratio 
show that indeed these liquid-crystalline phases can be 
stabilized solely from exclusion interactions, but no lay- 
ered ordering, of smectic (S) type, is found. However, the 
work of van der Kooij et al. [2| seemed to indicate that 
smectic or lamellar phases might exist in colloidal platelet 
suspensions at high volume fractions. Usually these sus- 
pensions are polydisperse in both thickness and diame- 
ter. The C phase stabilizes even for a considerable degree 
of diameter polydispersity, but uniformly thick platelets 
with very high diameter polydispersity cannot form a C 
phase, which presumably would be superseded by a S 



phase before the glassy or crystalline phases. Neverthe- 
less, the S phase remains elusive [2| and swollen lamellar 
phases has been observed to form, directly from a I sus- 
pension, only in a few materials [5|, [6| . 
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FIG. 1: (Color online). a-Zirconium phosphate platelets, (a) 
TEM of pristine ZrP platelets, (b) TEM of a ZrP monolayer 
platelet exfoliated by TBA. (c) Thickness schematic of a ZrP 
monolayer platelet. Length is not drawn to scale. The scale 
bars in (a) and (b) are 1 /Ltm. 

Here we report the observation of a stable S phase 
in a new system of platelets with high polydispersity in 
diameter but identical thickness. Contrary to the ap- 
proximately hard-core platelets analyzed in Ref. Q , our 
particles interact via long-range forces. We find that, as 
the platelet volume fraction <p increases, the samples fol- 
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low the phase sequence I-N-S, with no stable C phase 
in the volume-fraction interval explored. Some features 
of the phase transitions can be explained in terms of a 
density-functional theory (DFT) for a simple interaction 
model. 

Pristine a-Zirconium phosphate (ZrP) layered 
platelets were synthesized through a hydrothermal 
method similar to that reported in Ref. 0] and 
exfoliated with tetra-n-butylammonium hydroxide 
(TBA + OH~, Aldrich, 1 mol/1 in methanol) at a molar 
ratio of ZrP:TBA = 1:1 in de-ionized water. No salts 
were added during the experiments. After reaction, 
the platelets were allowed to expose in sonication for 5 
mins and then sit for one day to reach phase transition 
equilibrium. The platelets, which contain two TBA 
monolayers, have an identical thickness of L = 2.68 nm. 
This process has been demonstrated Q to lead to stable 
aqueous colloidal suspensions (see Fig. [T]). The average 
platelet is cylindrical but with relatively high variations 
in transverse shape. Samples containing various platelet 
volume fractions but the same polydispersity were 
prepared. Transmission Electron Microscopy (TEM) 
images [Fig. QJb)] allow for an approximate estimation 
of mean particle diameter (~ 2 (an). 
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FIG. 2: (Color online). Optical image (bottom) and images 
under cross polarizers (top) of samples. Bottom sequence 
corresponds to platelet volume fractions (from left to right) 
(j> = 0.004, 0.007, 0.009, 0.013, 0.017, 0.021, 0.031, 0.041, 
0.048, 0.054, 0.059, 0.074, 0.100, 0.160, and 0.200. Range of 
phase stability and phase coexistence are indicated by double 
headed arrows. 

To obtain a more accurate estimation of mean platelet 
diameter and polydispersity, Dynamic Light Scattering 
(DLS) was performed on the exfoliated ZrP monolayer 
platelets, yielding a platelet diameter of 2295 nm (close 
to the TEM estimation) and a polydispersity a D = 32%. 
To search for possible spatial ordering of the platelets, 
aqueous suspensions were directly characterized by X-ray 
diffraction (XRD). 

At 4> < 0.01, the platelet suspension forms a random 
stable dispersion in water (I phase). For samples with 
high <j) 1 birefringent features through crossed polarizers 
(Fig. [2]) reveal the occurrence of a transition to a nematic 
phase for <j) > 0.01. In the interval 0.01 < <j) < 0.03, I 
and N phases coexist. A full nematic phase is obtained 
for 0.03 < cj) < 0.06. The macroscopic phase separa- 




FIG. 3: (a) Nematic volume fraction (sample volume occu- 
pied by the N phase) as a function of platelet volume frac- 
tion. Symbols: experimental results. Lines: theoretical re- 
sults with <r D = 0.00 (dotted), 0.32 (continuous), and 0.52 
(dashed), (b) Theoretical platelet volume fraction as a func- 
tion of polydispersity for the I-N transition. Solid lines: I 
and N cloud curves. Dashed lines: I and N shadow curves. 
The coexistence region is shaded. 



tion is completed within 24 hours. Fig. [3] shows how 
the fraction of total sample volume occupied by the N 
phase behaves as a function of 4> (symbols). Nonlinearity 
in the two-phase region is due to fractionation effects: 
the diameter polydispersity of the platelets in the I and 
N phases are different (monodisperse suspensions would 
follow the lever rule and exhibit linear behavior). 

At high volume fraction of platelets, multiple scatter- 
ing peaks up to 5 orders in XRD experiments are ob- 
served (Fig. 0|) when </> > 0.06. These scattering peaks 
are the (001) family of layered arrangements of colloidal 
platelets (S phase). No XRD peaks were obtained for 
<\> < 0.06 (I or N phases). The S phase has been fre- 
quently observed in liquid crystals of rod-like molecules 
with layered structure. However, these organic smec- 
tic phases normally lack high-order XRD peaks due to 
their small domain size. For our colloidal platelets, since 
the size of the mesogen is of the order of hundreds of 
nanometers to microns, the domain size can easily reach 
tens of microns to millimeters. The presence of higher- 
order diffraction peaks and their relatively strong inten- 
sity (Fig. [3J is the result of the large domain size. Fig. 
[5] shows the optical microscopy image of the smectic liq- 
uid crystal observed in our system. A ribbon-like smectic 
mesophase around 50 yum wide and with a length of up to 
a few millimeters has been observed. Considering the size 
of the plate mesogens (around 2 /xm, Fig. [T]) , the optical 
microscopy image clearly demonstrates the assembly of 
our colloidal platelets into smectic layers parallel to the 
the plane of the image. 

Unfortunately it was not possible to directly measure 
the nematic, <pN, and smectic, 4>s, volume fractions (i.e. 
fractions of volume occupied by platelets in the respec- 
tive phase). An upper limit for <f>$ can be roughly es- 
timated by assuming smectic layers to consist of sin- 
gle platelet sheets in a random close-packed arrange- 
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FIG. 4: (Color online). XRD intensities as a function of scat- 
tering angle for various volume fractions. From top to bottom: 
A = 0.200, 0.160, 0.100, 0.074, 0.059, and 0.054. 




FIG. 5: Optical microscopy image of the smectic domains. 



ment, which gives <f>s — 4>2DL/d, 4>2D — 0.815 being 
the random close-packing area fraction of polydisperse 
hard disks with cr D = 0.3 @], and d the smectic period. 
In Table Q] values of (f>s are collected under the column 
labeled max. 

In the experiments the S phase occupies a small, slowly 
expanding, region of the sample. It seems that the phase 
transition (which may take many days to complete due to 
the high viscosity of the samples) is not proceeding fast 
enough to complete within the finite observation time. 
The weak diffraction peaks in the XRD spectra for the 
first smectic state lend support to a weakly first-order or 
continuous scenario for the N-S transition. 

One interesting feature is that the layer spacing (Ta- 
ble Q} exhibits a large variation (close to three-fold) in 
the volume-fraction interval explored, which may point 
to the existence of long-range platelet repulsion due 



TABLE I: Characteristics of the smectic phase measured in 
the experiments. A<f> = <j> - c/>* cxp , with c/>* cxp ~ 0.059 the 
estimated experimental volume fraction at the spinodal point. 
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d/L 


L e f(/L 
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estim. (max.) 


0.059 


0.000 


9.36 


3.49 


2.67 


0.108 (0.233) 


0.074 


0.015 


7.89 


2.94 


2.27 


0.136 (0.277) 


0.102 


0.043 


6.17 


2.30 


1.80 


0.193 (0.354) 


0.161 


0.102 


4.39 


1.64 


1.34 


0.326 (0.497) 


0.199 


0.140 


3.54 


1.32 


1.11 


0.441 (0.617) 



to platelet surface charges or possibly hydration forces 
(since platelets are hydrophilic, water-mediated effects 
might be at work 0]). Attractive van der Waals forces 
may also operate at shorter distances. 

To rationalize these findings, we map the real particle 
onto an effective model system of parallel polydisperse 
hard cylinders of diameter D and thickness L g- > L. 
Essentially, the model incorporates long-range repulsion 
via an effective particle thickness (a procedure used to 
model aqueous susp ensions of charged virus particles 
with added salt [10j]). For the moment we simply note 
that the value of L e ff will depend on <p and the type of 
phase. 

We obtained some predictions for the hard-core model 
by applying a theoretical analysis based on DFT. In the 
case of the I-N transition we used Onsager theory [i"ll |. 
which takes into account excluded-volume effects at sec- 
ond order in density and should give accurate results 
because of the high diameter-to-thickness ratio of the 
platelets. Since samples are polydisperse, a given dis- 
tribution of platelet diameters has to be imposed. Fur- 
thermore, the nematic order parameter, measuring the 
degree of orientation of the particle main axis about the 
director, is expected to be very high and the simplify- 
ing Zwanzig approximation jl2j (where only three mutu- 
ally perpendicular particle orientations, along x, y and z 
Cartesian axes, are permitted) can be implemented. The 
total free-energy density in reduced thermal units (fcsT), 
$ = ^id + $exc, can be split into ideal, 

$ id = J dr{2 P± (r) pnp x (r)-l]+/>||(r) [Inpy (r) - l] } , 

(1) 

and excess 
3>cxc(7r(i}) 3 /4) = 2 



(1) (2) , (1) (2) , (2) (1) 

rrij_ m\ + m,, mj_ + ml m\ 



(2) 

parts. The moments of the number density distri- 
bution functions p M (r) are defined as 

m[f) - (n(D) 3 /4) J drr a p^r). (3) 

We use the dimensionless platelet diameter r — D/(D) 
(with (D) the mean diameter) as polydispersity vari- 
able, and used the notation p x (r) = p y (r) = p±(r) and 
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Pz(r) — P\\(r) to label the densities of species oriented 
perpendicular and parallel to the director, respectively. 
For the I and N symmetries we have 



Px(r) = p\(r) = 



p\r) 



(4) 



Pl(r) = 1 —§^-p N (r), ^(r) = I±^p», 

(5) 

where —1/2 < Q(r) < 1 is the N order parameter of 
species with reduced diameter r. To find the I and N 
coexistence densities, i.e. 



„i,n 



fi—x,y,z 



drp l ^(r), 



(6) 



in a system with a relative fraction of volume occupied by 
the N phase given by the parameter < 7 < 1, we need 
to implement the constrained functional minimization of 



$ = (1 - 7 )$ I +7$ f 



(7) 



with respect to pJ; N (f*). are the total free-energy 

densities of the coexisting I and N phases. The constraint 
is dictated by the lever rule restriction: 



p (r)^(l- 1 )p 1 (r)+ 1 p N (r), 



(8) 



where po(r) — poh(r) is the parent number density dis- 
tribution function with a fixed function h(r), which is 
the diameter distribution function density normalized 
as J drh(r) = 1. We choose a two-parameter function 
which, after normalising and applying the additional con- 
dition J drrh(r) = 1 (note that this is equivalent to 
(r) = 1), becomes 



h(r) = 



2T[(i 



r[(z, + i)/2f 



2)/2]" +1 „ 
2~r ex P ' 



r[fr + 2)/2] 
r[(i/ + i)/2]' 



(9) 

where T(x) is the Gamma function. The free parameter 
v > controls the degree of polydispersity. Note that 
the present function h(r) has a Gaussian tail for r ^> 1, 
while it fulfills h(Q) = 0, in contrast to the Gaussian 
distribution function. The constrained minimization of 
Eq. ([7]) and the equality of osmotic pressures 



I,N 



P 



I.N 



I.N 



(10) 



between both phases allow us to find the coexistence total 
densities as a function of 7. For 7 = we find the coex- 
istence between a I phase which occupies the entire vol- 
ume, called the cloud phase, and an infinitesimal amount 
of the N phase (called shadow phase). In the other limit, 
7 = 1, we find the cloud N-shadow I coexistence. Results 
for the I-N transition are plotted in Fig. EJa) as lines for 
three different polydispersities: er D = 0.00 (monodisperse 
fluid), 0.32 (same as in experiment), and 0.52. The ef- 
fective thickness (assumed to be independent of <p) was 
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FIG. 6: (a) Theoretical volume fraction of the spinodal in- 
stability against smectic fluctuations as a function of polydis- 
persity. (b) Osmotic pressure (in reduced units) vs. platelet 
volume fraction cj> for smectic (continuous line) and nematic 
(dashed line) for <j d = 0.52. Filled circles: smectic states; 
open circles: unstable nematic states with <j> obtained as 
(j> = (j>* thc + A(j>, with <^* tho = 0.288 (theoretical bifurcation 
point) and A<f> taken from experiment (see Table |TJ. (c) Den- 
sity profile p(z,r) for r = 1.5 (solid line), 0.4 (dashed) and 0.8 
(dotted), (d) Normalized density distribution function h(z,r) 
in one smectic period. Lines as in panel (c). In (c) and (d) 
the smectic volume fraction is <j>s ~ 0.452 and the smectic 
period d/L c g — 1.211. 



set to -Leff = 9L so as to make the experimental and the- 
oretical curves with the same value of er D approximately 
overlap (note however that our theoretical data, with- 
out scaling, agree well with the experimental results of 
van der Kooij et al. for their approximately hard 
platelets). We note that: (i) for high cr D the curvature in 
the transition region is similar to that of the experiment, 
(ii) The I+N coexistence gap is too small in the theory; 
this may be due to the crudeness of the model, and/or to 
an inaccurate experimental estimation of er D (note that 
the width of the I+N region in the theoretical curve for 
<7 D = 0.52 is quite close to the experimental one). As in- 
ferred from Fig. Eljb), the coexistence gap (defined as the 
distance between the I and N cloud curves at the same 
polydispersity) increases with a u due to fractionation ef- 
fects: larger platelets tend to go to the N phase, while 
smaller platelets preferentially populate the I phase. 

Next we turn to the theoretical description of the N-S 
transition. The S phase demands a description in terms 
of the density distribution p(z,r), where z is the dis- 
tance along the nematic director. From this, the lo- 
cal smectic volume fraction can be calculated, 4>s{z) — 
vq J drr 2 p(z, r) (with vq = it (D) L e g/4 an average par- 
ticle volume), and also a mean smectic volume fraction, 

4>s = d Jq dzcf>s(z). A DFT fundamental-measure 
scheme for strictly parallel hard cylinders [131 ] was gener- 
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alized to the polydisperse case (the assumption of parallel 
particles should be valid since the nematic order parame- 
ter is expected to be very high). The excess, z-dependent 
local free-energy density in units of fc^T is 



^oxc «o = -n ln(l - n 3 ) + 



l-«3 {l-n 3 f 
(11) 



where the weighted densities are 

n (z)=M (0) *w(°)(z), n 3 (z) = X (2) * w (3) (z), 
n^z) = ,M« * w<°)(z), nf (z) = M (0) * w (3) (^), 
n^(z) = M (1) *cj( 3 >(z), njj(z) = .M (2) * u/W(z), 

(12) 

with M {a) {z) = m ( f \z)/K and k = {D)/L cH . The 
weights are 



^ {Z ) = 1 S ( 1 --\ Z \), w o)(z) = e(i-w 



(13) 



with <5(z) and G(z) the usual Dirac-delta and Heaviside- 
step functions; the symbol * denotes convolution, i.e. 



1 



f*wW[z) = -[/(*- 1/2) + /(*+ 1/2)], 

rZ+1/2 

f* U &\z) = / dz'f(z'). 

J 2—1/2 



(14) 



In all preeceding equations (fT2|) - (fT4l) the coordinate z is 
in units of L e ff- The ideal and the interaction free-energy 
functionals per unit volume for the smectic are then 



d' 1 I dz 
la 



drp(z, r) [In p{z, r) — 1] . 



Pfexc/V ^d^ dz<t> exc (2 

Jo 



(15) 
(16) 



The equation of state of S phase was obtained by func- 
tional minimization of the total free-energy with re- 
spect to p{z,r) = pu(z,r) with the restriction po(r) — 

d~ l Jo dzp(z, r), po(r) being the (fixed) mean density 
distribution function per period. The N-S transition is 
found to be continuous for all values of polydispersity 
(see Fig. [B] (a), where the limit <r D = corresponds to 
the one-component model of Ref. [IH). This is in qual- 
itative agreement with experiments which, as discussed 
above, suggest a continuous or weakly first-order phase 
transition [X 51 ] - 

Fig. [IJb) shows the osmotic pressure for N and S 
phases as a function of tj> for the hard-core potential. 
The N phase is unstable beyond the bifurcation point. 
Even though there are no signs of N+S coexistence or 
macroscopic fractionation at the transition, some subtle 



micro-fractionation effects in the S phase exist, shown 
in Figs. (Hfc-d) which pertain to values cr D = 0.52 
and 4>s = 0.452. Density profiles for different diame- 
ters show micro-fractionation: the probability density 
h(z, r) — p(z, r)f J drp(z, r) shows that the relative frac- 
tion of large platelets is higher at smectic layers, while 
the opposite occurs for smaller platelets, with maxima 
at interstitials. The relative fraction of intermediately- 
sized platelets has two maxima per period, located a bit 
off the layers. 

Finally we extract some information from the model in 
order to interpret experimental data and, in particular, 
to refine the values for 4>s previously derived as upper 
estimates. We assume the experimental transition to oc- 
cur at a platelet volume fraction 4>* m 0.059 (interpreted 
as a spinodal point or as the point where the smectic 
phase appears for the first time). The experimentally 
observed N states for <fi > <fi* would then be unstable, 
with nematic volume fractions 4>n < 4>s'i here (f>s is the 
(unknown) platelet volume fraction of the stable smectic 
state (at the same osmotic pressure). Since the smectic 
volume in the samples is small, we take c/)n — 4>- Also, 
we identify the (unstable) experimental nematic states 
with the theoretical (unstable) nematic branch. We pre- 
dict <f>s from our equation of state, Fig. [|3b), as follows. 
First, the volume-fraction distance A<fi = <f> — <f>* is taken 
to be the same in both theory and experiment, which 
allows to obtain, for each <p (or Acj)) in Table HI a theoret- 
ical value [indicated by filled circles in Fig. EUb)], to 
which our hard-core theory uniquely associates a value 
d/L e g, and estimate L c r/L = (d/L)/(d/L e g) (see Table 
H|); from low to high volume fraction, L e g/L — 2.67-1.11. 
Now the estimated (f>s for the real system follows from 
4>s = 0g he /(Leff/i) -see Table U under column estim. 
[161] . Even though the refined procedure probably gives 
more reasonable values for (j>s than the previous crude 
estimates, we cannot at present improve the model for 
lack of accurate experimental measurements. However, 
all evidence points to the existence of a weak phase tran- 
sition from a N to a stable S phase in this novel colloidal 
suspension of platelets. 

To summarize, we have experimentally analyzed a new 
highly polydisperse colloidal platelet suspension exhibit- 
ing smectic ordering, thus confirming that this elusive 
phase can be stabilized in platelet colloids when poly- 
dispersity suppresses the columnar phase. A theoretical 
model for polydisperse platelets based on DFT and an 
effective hard-core model partially explains the experi- 
mentally observed phase behavior. 
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